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Abstract. We study the Brownian force model (BFM), a solvable model of avalanche statistics for 
an interface, in a general discrete setting. The BFM describes the overdamped motion of elastically 
coupled particles driven by a parabolic well in independent Brownian force landscapes. Avalanches 
are defined as the collective jump of the particles in response to an arbitrary monotonous change 
in the well position (i.e. in the applied force). We derive an exact formula for the joint probability 
distribution of these jumps. From it we obtain the joint density of local avalanche sizes for stationary 
driving in the quasi-static limit near the depinning threshold. A saddle-point analysis predicts the 
spatial shape of avalanches in the limit of large aspect ratios for the continuum version of the 
model. We then study fluctuations around this saddle point, and obtain the leading corrections 
to the mean shape, the fluctuations around the mean shape and the shape asymmetry, for finite 
aspect ratios. Our results are finally confronted to numerical simulations. 
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1. Introduction 

A large number of phenomena, as diverse as the motion of domain walls in soft magnets, fluid 
contact lines on rough surfaces, or strike-slip faults in geophysics, have been described by the 
model of an elastic interface in a disordered medium HI Si. A prominent feature of these 
systems is that their response to external driving is not smooth, but proceeds discontinuously by 
jumps called “avalanches”. As a consequence of this ubiquitousness, much effort has been devoted 
to the study of avalanches, both from a theoretical and an experimental point of view in 0, El E] . 
Despite this activity, there are few exact results for realistic models of elastic interfaces in random 
media. 

An exactly solvable model for a single degree of freedom, representing the center of mass of 
an interface, was proposed by Alessandro, Beatrice, Bertotti and Montorsi (ABBM) [8l E] on a 
phenomenological basis in the context of magnetic noise experiments. It describes a particle driven 
in a Brownian random force landscape. In mm it was shown that for an elastic interface with 
inhnite-ranged elastic couplings, the motion of the center of mass has the same statistics as the 
ABBM model. 

In this article, we study a multidimensional generalization of the ABBM model, the Brownian 
force model (BFM). This model, introduced in [TTl fl^ [131 E]; was shown to provide the correct 
mean-held theory describing the full space-time statistics of the velocity in a single avalanche for 
d-dimensional realistic interfaces close to the depinning transition. Remarkably, restricted to the 
dynamic of the center of mass, it reproduces the ABBM model. This mean-held description is 
valid for an interface for d > due with due = 4 for short ranged elasticity and due = 2 for long 
ranged elasticity. 

As shown in mm the BFM has an exact “solvability property” in any dimension d. It is 
thus a particularly interesting model to describe avalanche statistics, even beyond its mean-held 
applicability, i.e. for any dimension d and for arbitrary (monotonous) driving. It allows to calculate 
the statistics of the spatial structure of avalanches, properties that the oversimplihed ABBM model 
cannot capture. In Ref. na some hnite wave-vector observables were calculated, demonstrating 
an asymetry in the temporal shape. Very recently the distribution of extension of an avalanche 
has also been calculated na. 

In this article we study a general discrete version of the BFM model, i.e. N points coupled 
by an elasticity matrix in a random medium, as well as its continuum limit. In the discrete model 
each point experiences jumps Si upon driving. We derive an exact formula for the joint probability 
distribution function (PDF) T’[{5'j}] of the jumps Si (the local avalanche sizes) for an arbitrary 
elasticity matrix. In the limit of small driving this yields a formula for the joint density p[{5'j}] 
of local sizes for quasi-static stationary driving near the depinning threshold. This allows us to 
discuss the “inhnite divisibility property” of the BFM avalanche process. The obtained results 
are rather general and contain the full statistics of the spatial structure of avalanches. They 
are, however, difficult to analyze in general since they contain many variables, and thus require 
computing marginals (i.e. probabilities where one has integrated over most of the variables) from a 
joint distribution. This is accomplished here in detail for the fully-connected model. We hnd that 
in the limit of large N there exist two interesting regimes. The hrst one corresponds to the usual 
picture from mean-held depinning models mm, whereas the second one is novel and highlights 



Avalanches in the BFM 


3 


the intermittent nature of the avalanche motion. 

We then analyze the shape of avalanches, hrst in a discrete setting by considering few degrees 
of freedom. The probability exhibits an interesting saddle-point structure in phase space. We then 
study the continuum limit of the model. We find that the spatial shape of avalanches of fixed total 
size S and extension i, becomes, in the limit of a large aspect ratio S/i"^, dominated by a saddle 
point. As a result, the avalanche shape becomes deterministic, up to small fluctuations, which 
vanish in that limit. We calculate the optimal shape of these avalanches. We then analyze the 
fluctuations around the saddle point. This allows us not only to quantify the shape fluctuations 
seen in numerical experiments, but also to obtain the mean shape for avalanches with smaller 
aspect ratios. We test our results with large-scale numerical simulations. While our results are 
obtained in the special case of an elastic line with local elasticity {d = 1 ) the method can be 
extended to other dimensions d and more general elasticity. Finally, we discuss the applicability of 
our results to avalanches in realistic, short-ranged correlated disorder. The outline of this article 
is as follows: Section [^recalls the definition of the BFM model, which is first studied in a discrete 
setting with general, non-st at ionary driving. The results of [121 HSl IE] allow us to obtain the 
Laplace transform of the PDF of local avalanches sizes. Section contains the derivation of the 
main result: the full probability distribution of the local avalanche sizes. Section focuses on 
the limit of small driving, and how to obtain the avalanche density. Section contains a detailed 
analysis of the fully-connected model. Section studies avalanche shapes for interfaces with a 
few degrees of freedom. Section contains one important application of our result, namely the 
deterministic shape of avalanches with large aspect ratio for an elastic line. Section analyses 
the fluctuations around this optimal shape. Section discusses the application of our results to 
short-ranged disorder and quasi-static driving. A series of appendices contains details, numerical 
verifications and some adjunct results. In particular, in Appendix C we introduce an alternative 
method, based on backward Kolmogorov techniques, to calculate the joint local avalanche-size 
distribution, following a kick in the driving. 


2. The Brownian force model 

2.1. Model 

We study the over-damped equation of motion in continuous time t of an “interface”, consisting 
of N points with positions G M, z = 1,..., A^. Each point feels a static random force Fi{uit) 
and is elastically coupled to the other points by a time-independent symmetric elasticity matrix 
Cij with Dj = 0. Each particle is driven by an elastic spring of curvature m? centered at the 
time-dependent position Wa. The equation of motion reads 

N 

vdtUit = ^ CijUjt - - Wit) + Fi{uit) ( 1 ) 

i=i 

for z = 1... A^. The Fi{u) are N independent Brownian motions (BM) with correlations 


[Fiiu) — Fi{u')Y = 2a\u — u'\ , Fi{u)Fj{u') = 0 for z 7 ^ j 


( 2 ) 
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and Fi{u) = 0; the overline denotes the average over the random forces Fi{u). For dehniteness we 
consider [|] a set of one-sided BMs with u>0 and Fi{0) = 0. 

We furthermore suppose that (i) the driving is always non-negative: Wt,i, wu > 0, and (ii) 
the elastic energy is convex i.e. c^- > 0 for i ^ j. Under these assumptions, the Middleton theorem 
[TB] guarantees that if all velocities are non-negative at some initial time: 3fo ^ IR|Vi, Uhq > 0, 
they remain so for all times: > to, iiu > 0. 


Some explicit examples of elasticity matrices: Throughout the rest of this article, we sometimes 
specify the elasticity matrix. The models studied are (where c denotes the elastic coefficient): 

(i) The fully connected model: = c(U — Stj) 

(ii) The elastic line with short-range (SR) elasticity and periodic boundary contitions (PBCs) 
Cij = c {Sij-i + 6i-ij - 2dij) with i + N = i 

(iii) The elastic line with SR elasticity and free boundary conditions: 

Cij c [(5j -|- 6i—i j 6ij(f2 5i\ <^i7v)] 

(iv) The general d-dimensional elastic interface with PBCs, where i and = c(/(||i — i||) — 

~ill))j here ||i — j|| is the Euclidean distance in lA and /(r) the elastic kernel. 
Long-ranged elasticity (LR) is usually described by kernels such that /(r) ~ ^ 

in Fourier). 


2.2. Velocity Theory 


Supposing that we start at rest for f = 0, Uij=o = 'Ui,t=o = 0, then it is more convenient (and 
equivalent) to study the evolution of the velocity held directly. The equation of motion reads 

N 

pdtiiit = ^ Cijiijt - - wu) -f \/2aUitCt > (3) 

i=i 


where the Q are N independent Gaussian white noises, with = 5{t—t')5ij and = 0. Equation 
([^ is taken in the ltd sense. Note that we replaced the original quenched noise dtFiiuu) by an 
annealed one \/2aUitQ, making Eq. ([^ a closed equation for the velocity of the interface. The fact 
that ([^ and (|^ are equivalent (in the sense that disorder averaged observables are the same) is a 
non-trivial exact property of the BFM model. It was hrst noted for the ABBM model nu and 
extended to the BFM [T^ ITT] . It originates from the time-change property of the Brownian motion 
dB{f{t)) =in law \/ f'{t)dB{t) for increasing f{t) = ut, valid as a consequence of the Middleton 
property Ut>0. A derivation of this property is recalled in [Appendix A 


2.3. Avalanche-size observables 

In this article we focus on the calculation of avalanche-size observables dehned in the following 
way. Starting from rest at f = 0 as previously described, we apply a driving wu > 0 for f > 0 
during a hnite time interval such that f^dtwu = Wi (stopped driving protocol). In response to 

f The model can also be studied in a stationary setting, see e.g. [13 [H]. 
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this driving, the points move and we dehne the local avalanche size Si as Si = dtitu, that is 
the total displacement of each point. We adopt the vector notation 

S = {Si,...,Sn) , w = {wi,... ,wn) ■ (4) 

The SiA are random variables whose statistics is encoded in the Laplace transform, also called 
generating fnnction G'(A), and dehned as 

G{X) ='^ . (5) 

The BFM possesses a remarkable “solvability property” that allows ns to express this fnnctional 

as [I3l[ll] 

G(A) = = e™' ^-1 (6) 

in terms of the solntion Ui of the “instanton” eqnation. The latter reads 

N 


Xi = 


-aul + CijUj , 

i=i 


where we have defined the dimensionless matrix 

1 


Cij 6ij rjCij 


(7) 


( 8 ) 


which contains all elastic and massive terms in the instanton equation. The solution of Eq. Q 
which enters into Eq. ([^ is the unique set of variables Ui continuous in Xj with the condition that 
all Ui = 0 when all Xj = 0. The derivation of this property is recalled in a discrete setting in 


Appendix A The instanton equation thus allows us in principle to express the PDF P{S) of the 


local avalanche sizes, as the inverse Laplace transform of G(A). In the next section we obtain P{S) 
directly, without solving Q, which admits no obvious closed-form solution. We will note (...) the 
average of a quantity with respect to the probability P. Note that the PDF P{S) depends only 
on the total driving Wi = dtwit and not on the detailed time-dependence of the Wu- This is a 
particularity of the BFM model. 


2.4- The ABBM model 

Before going further into the calculation, let us recall the result of Ref. that the statistical 

properties of the center of mass of the discrete BFM model is equivalent to that of the ABBM 
model. To be precise, if we write the total displacement (i.e. swept area) u* = Yhi'^n total 
drive = 'Y4i Wu then, in law, we have 

pdtiit = - vJt) + \/2aut^t ■ ( 9 ) 

Here is a Gaussian white noise ^t^t' = d{t — t') and = 0. [§] This equivalence implies that the 
PDF of the total avalanche size S = dt Ut = i^ discrete BFM model, following an 

arbitrary stopped driving dt w* = w, is given by the avalanche-size PDF of the ABBM model 

PEI [13], 

§ Note that this result uses = 0 that the center of mass obeys the same equation with a noise scaled as 

and driving by 
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Here Sm is the large-scale cutoff for avalanche sizes induced by the mass term. This first result 
on a marginal of the joint distribution P{S) will provide a useful check of our general formula 
obtained below for N > 1. 

3. Derivation of the avalanche-size distribntion in the BFM 


For simplicity we now switch to dimensionless units. We define 


a 


Wi 






Sr, 


Q - A 

— Q 5 

^m. 


where Sm = The instanton equation ([^ now reads 

N 

Xi = —vl + ^ Cr 


yrjVj . 


( 11 ) 


( 12 ) 


i=i 


The generating functional is given by 


G(A) = G(A) = . (13) 

In the following we drop the tildes on dimensionless quantities to lighten notations, and explicitly 
indicate when we restore units. For the ABBM model, it was possible to explicitly solve the 
instanton equation for the generating function G{X). The inverse Laplace transform was then 


computed, leading to (10). Here this route is hopeless because Eq. (12) admits no simple closed- 


form solution. We instead compute directly the probability distribution P{S) using a change of 
variables in the inverse Laplace transform (ILT): 

N „ 

/ d^Aexp (-A 


= I ^ 


■ G'(A) 


(14) 


1 

2i7r 


N 


= / dm 


duAT det 1 7 -T ) exp 
\dvj' 


N 


E( 


N N 

CijVj^Si ^ ^ VjWi 
i=l j=l i=l 

where “i” denotes the imaginary unit number to avoid confusion with indexes. The first formula 
is the ILT where we left unspecified the multi-dimensional contour of integration C. In the second 
line we used the expression of Aj in terms of Vj from (12), as well as the dimensionless version of 


(§. Changing variables from Aj to Vj, the contours of integration are chosen to obtain a convergent 


integral, see second line of Eq. (14). This makes this derivation an educated guess, which however 


det|'A| = 

dvj 


is verified in Appendix B We also give another derivation for a special case in Appendix C To 
pursue the derivation, the Jacobian is written using Grassmann variables as 

. N / N \ 

JJ diJidipi exp I ^ j . (15) 

i=l V*J=1 / 

Reorganizing the order of integrations and changing m —)■ im, we write 

= ( X- ) n / d^id'i/’i Yl / dm exp ( - ^(n^ + ^ iCijVj)Si 

i=iJ V i=i 3=1 

N N 

+ ^ iViWi + ^ 

*P=i 


2=1 


(16) 
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Integrating on Vi leads to 
' 1 ^ '■ 


N 


N 


2ti J 


n / d*d*(i)(''/=i ms* “p -7E 


2'ijjitjji - CijSjY 


i=l 


v2=l 


Si 


(17) 


2=1 
N 

+ ^ ^ 'tpiCij'tpj 
7J = 1 

Finally, using 'ipf = 'll)j = 0, the integration over the Grassmann variables can be expressed as a 
determinant, leading to our main result 

.N / N 


P{S) = 


2^ 


1 JL [Wi 


nSi exp 


n c 'iS' 

1=1 


. 2=1 


2=1 


5,: 


det (Mjj) 


NxN 


(18) 


- Ylk=l ^ikSi 

' s. 


Mij — Cij + Si 


Cij Sij Cij 


Here c^- is the elasticity matrix. This is the joint distribution expressed in dimensionless units 
(0. The expression in the original units is recovered by substituting Si —)■ Si/Sm, Wi —t Wi/Sm 
and P —>■ S^P in (18) while keeping Cij hxedfjj] 

Note that for zero coupling, = 0, Eq. (18) becomes P{S) = Hili -PA.BBM(*S'j): the different 
points are decoupled and one retrieves N independent ABBM models. Non-trivial tests of the 
formula are performed in Appendix B One general property is that the average local size is 
i^i) — This average gives the shape of the interface in the large-driving limit. When 

Wj S> 1 uniformly in i, it is easy to see by expansion of the above formula that Si = {Si)+0{^/7Fi)rji 
where r]i are (correlated) Gaussian random variables. 

We show in [Appendix G using different methods, that when the driving is in the form of 
kicks. Wit = WiS{t) P{S) satisfies the exact equation 


N 

E 

q :=1 


dP d^P dP 

} CajWj + - Wn 


dw, 


i=i 


dwl, 


‘dS^ 


= 0 . 


(19) 


We also show that ( |18| ) solves this equation. This alternative derivation support our result (18) 
ans shed some light on its structure. 

Interpretation: Some features of our main result can be understood as follows. Gonsider the 
equation of motion ([^. Upon integration from t = 0 to t = oo we obtain 

^ roc 

0 = ^ CijSj - m‘^{Si -Wi)+ / dt^/2^tCt ■ ( 20 ) 

i=i 


If we could replace the sum of white noises by a gaussian random variable 

poo / poo _ 

/ dt <^2aUitCt -^ \ 2a / dtuu E, ^ p2aSi E, , 


( 21 ) 


II Note that this formula can be generalized to the case of site-dependent masses and disorder strengths, mi^ar. 
the expression in the original units is obtained by the substitution Si —)■ S'i/S'm, Wi —)■ Wi/S^^ and P —)■ ]([j Sl^P in 


(18|) with Sl„ = ^ and Cu = 5^ -U cn- 

'* V m rn- ‘'J ‘'J m- ‘'J 


*rThis is sufficient, since we noted above that the result does not depend on the detailed time-dependence of the 
driving. 
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then we would obtain ( |l8| ), but with a slightly different determinant given by the replacement 
5ij —>■ ^5ij in Mij in (h^. However, the replacement (21) is not legitimate because the variables 


iiit are correlated in time. The determinant in (18) takes care of that correlation. 


Probability distribution of the shape 
Section 


2.4 


Even if it is far from being obvious on Eq. (|18|), we know from 


that the probability distribution of S' = Xlili given by (10) with w = 


Wi 


This allows us to define the probability distribution of the shape of an avalanche, given its total 
size S: Consider si,..., s^v G [0,1] with sat = 1 — such that Si = Ssi. The probability 

distribution of the s* variables, given that the avalanche has a total size S = 

\2N ^ 


CW+I 

P(s|S') = 2\/Tr -exp 


w 


(S — w)" 

4^ 


P(Si) 


2 = 1 


Si 1 . 


( 22 ) 


4. Avalanche densities and quasi-static limit 

The goal of this section is to dehne and calculate avalanche densities. These allow us to describe 
the intermittent motion of the interface in the regime of small driving, Wi small. The dependence 
of the PDF, P^{S), on the driving is denoted by a subscript w. We first study the jumps of the 
center of mass described by the ABBM model. 


f.l. Center of mass: ABBM 


For the ABBM model (and for the total size S = f^e BFM model) the avalanche-size 

PDF is given by 


Pw(^) 


W 

y exp 
2 





(23) 


where w = total driving. The limit of small driving w is very non-uniform. In the 

sense of distributions, its limit is a delta distribution at S' = 0, 


P^{S) 5{S) . (24) 

However, this hides a richer picture and a separation of scales between typical small avalanches 
S ~ and rare large ones S' ~ 1. If one defines S' = w^s, the PDF of s has a well-defined w —> 0 
limit given by 

Po{s) = ^ 3 exp () , (25) 

2^/tts2 V 

which is indeed normalized to unity / dspo{s) = 1. Hence avalanches of sizes S ~ are typical 
ones. However, all positive integer moments of po are inhnite. This indicates that these small 
avalanches, though typical, do not contribute to the moments of which are hnite and controlled 
by rare but much larger avalanches which we now analyze. In the limit of small w, there remains 
a probability of order w to observe an avalanche of order 1. For fixed S' = 0(1) ^ one has 

Pw(P) = wp(P) + 0(w2) , p{S) = exp (^-0 . (26) 
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This defines the density (per nnit w) of avalanches. These are the “main” avalanches with S' S> w^, 
which are also called “qnasi-static” avalanches (see below and Section [^. The density is not 
normalizable becanse of the divergence at small S, bnt all its integer moments are finite and 
contain all the weight in that limit, i.e. (S"") = w JdSp(S)S"' + O(w^). In particnlar, (S) = w 
implies f dSp(S)S = 1. 

We now show that the avalanche density contains more information and controls the moments 
even for finite w, a property that follows as a consequence of Pw{S) being the PDF of an infinitely 
divisible process. This is best seen on its Laplace transform 

Gw(A) = J dSe^^P^iS) = , Z{X) = ^(1 - ^/^^) . (27) 

The “infinite-divisibility property” indeed follows: Vm and Vw = wi Wm such that Wi > 0 

m 

Gw(A) = nG'wi(A) , Pw(5) = (Pwi*---*PwJ(5) , (28) 

i=l 

where * denotes the convolution operation. Hence S' is a sum of m independent random variables 
for all m. The ABBM avalanche process can thus be interpreted as a Poisson-type jump 
process (a Levy process) with jump density p{S) [1^. In general the density can be defined 
as p{S) = |w=o for fixed S' > 0 (i.e. it does not hold in the sense of distributions), and the 

relation between Z{X) := |w=o and p is 

Z(A) = j d^(e^^ - l)p(S) . (29) 


The —1 takes care of the divergence at small S'. This allows us to write the relation between P^ 
and p, expanding (27) in powers of w, as 

r Yi f* 

I dSe^'^Pw(S') = X] ^ / dsi ■ ■ ■ ds„(e^^^ - 1) ■ ■ ■ (e^^" - l)p(si) ■ ■ ■ p(s„) 


n=0 


(30) 


Taking derivatives w.r.t. A, this decomposition shows that the (positive integer) moments of Pw are 
entirely controlled by p, for arbitrary fixed w (beyond the small-w limit). In this sum the term of 
order w”' can be interpreted as the contribution to the total displacement S' of the interface (after 
a total driving w) of a n-avalanche (quasi-static avalanche) event (of order 0(1)). The convolution 
structure in (30) shows that these events are statistically independent in the ABBM model. In 
this model however, this interpretation only holds at the level of moments. The accumulation of 
infinitesimal jumps, manifest in the non-normalizable divergence of p at small S' prevents us to 
extend this interpretation to the probability itself, see [Appendix D| for a discussion. 


4 . 2 . BFM 

In the BFM, “the infinite-divisibility property” of the avalanche process is even richer, since 
avalanches occur at different positions along the interface. Let us define the j-th “elementary” 
driving which applies only to site j, i.e. Wi = Wj 6 ij, and denote the corresponding size-PDF as 
Pu,^.(S'). Consider now the PDF for the general driving, P^(S'). From the structure of its LT, see 
(13), as a product of exponential factors linear in the Wi, this PDF can be written as a convolution 
for w = {wi, ...,wn), 

P^S) = Pu„{S) * ■ ■ ■ * P^^{S) . 


(31) 
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An avalanche in the BFM can thus be understood as a superposition of N avalanches independently 
generated by each local driving Wj. 

As for the ABBM model (center of mass), the structure of the LT of the PDF P^^(S') shows 
that each of these elementary jump processes is inhnitely divisible. We dehne the avalanche density 
generated by the driving on the j-th point as 


^ _ dP^iS) dP^^iS) 

Pj[o) ■— - — -1,^=0 


dwj dwj 


|ui,=0 


where as in the previous case, this equality is to 
Consider the functions Vj of A which appear in Eq. 


(32) 

De understood point-wise in the S variables. 
13) and satisfy Eq. (12). It is the analogue of 


Z{X) appearing in (27) for the ABBM model and we t 

V. = f d^S 


ins conjecture the generalization of (29), 




1 ft(S) . 


(33) 


This allows us to write an equation relating Pwj{S) to Pj{S) similar to (30) (see Appendix D). 


The subtleties linked with the accumulation of small avalanches and the non-normalizability of 
Pj{S), are the same as in the previous case, which is also reminiscent of the fact that the limit 
of small driving of Pw{S) is very non-uniform, as we now detail. Consider Wi = wfi with w —)• 0 
and fi hxed: the limit of P^{S) is again given (in the sense of distributions) by n^i^(‘5'i)- More 
precisely, in this small-tc regime, almost all avalanches are O(w^): Si = wfsi with the s* distributed 
according to 

N 

Po(s) = JJpo(si) , (34) 


i=l 


as can be seen from an examination of (18) in that regime. The PDF po was dehned in (25). One 
sees that the regime Si ~ contains all the probability, and that for these very small avalanches 
the local sizes are statistically independent. 

The remaining 0{w) probability to observe large avalanches Si = 0(1) is encoded in the 
densities Pj{S), 


N 


PdS) = y^w,p,iS) + Oiw^) 


i=i 


(35) 


As before, the positive integer moments are entirely controlled by pj. A more general expression, 
which illustrates that these large avalanches occur according to a Poisson process, is given in 
Appendix D 

We now give exact expressions for these densities. For a general elasticity matrix, the 


expression of pj is obtained from Eq. (18), and contains a determinant. Remarkably, one can 
compute this determinant in various cases, leading to the following result 


- (. 2 ^) 


N 


s. 


N 


(nf=i ^0 


T^(q\ I 1 (^7 = 1 C'p'S'j)" 

jK{S)e-x.Y> \-jy- 


2=1 


s, 


( 36 ) 


where K{S) depends on the chosen elasticity matrix: 


Fully connected model: K{S) = 


7V-1(EAi^ 


nf=i s. 
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Linear chain with periodic boundary conditions: K(S) = ^ 5 . 5 .^^ 

Linear chain with free boundary conditions: K{S) = ^ 


SiSn 


PDF of the shape in the small-driving limit As we just detailed, the small-driving limit of Pw{S) 
exhibits a complicated structure due to the accumulation of small avalanches. The situation is 
very different for the PDF of the shape of the interface conditioned to a given total size S = 0(1) 


(22). This conditioning naturally introduces a small-scale cutoff that simplihes the small driving 


limit Wi = wfi with tc —)• 0 which reads 

cV- 1-5 

p{^S) = limP(^5) = 20f^^exp 

^^0 l^ifi V 4 


5 


N 


i=i 


(37) 


This limit holds in the sense of distributions, and p(^*S') dehnes a normalized probability 
distribution. This indicates that the only small-scale divergence present in pj originates from 


the direction S'j ~ S' —)■ 0 uniformly in j, in agreement with the conjecture (33). 


5. Fully-connected model 


In this section we use our result (18) and analyze it for the fully-connected model with uniform 


driving. Most calculations are reported in Appendix E where we also consider driving on a single 
site, Wi = Wi6ii. 


Structure of the PDF and marginals In the fully-connected model with homogeneous driving 


Wi = w, it is shown in|Appendix E|that our main result ([1^ has the simple structure 
P(S) = 


w 


N 


w - 1 - cS/N 


]^Pto,S/w(*S'i) . 


i=l 


We dehned 


Pw,z{,Si) 


w -\- cz 


exp 


{w + cz - {IF c)Sj)" 
4P 


(38) 


(39) 


For each w, z > 0, it is a probability distribution, that corresponds to the (dimensionless, with 
= 1) PDF of the avalanches of one particle in a Brownian force landscape (ABBM model), 
interacting with one parabolic well through the force m^{w — Ui) and with another parabolic well 


through the force c{z — uf). Formula (38) is thus reminiscent of the fact that the various sites 


interact with one another only through the center of mass of the interface. This simple structure 


permits a direct evaluation of various marginals of (38) of the type P{{Si ,..., S'p}, S) (local sizes 


on p < A^ sites and total size). This is done in Appendix E Here we focus on the joint PDF of 
the total size S, and the single-site local avalanche size Si < S. Its explicit form is 


P{Si,S) = 


w 


2y/^Sf 


X exp 


-(iV-1) 


w + cS/N 

2V^{S-Si)^P 


exp 


{w + cS/N - (1 + c)^i)' 


45i 


(40) 


((AT -l){w + cS/N) - (1 + c)(^ - ^ 1 ))' 
4(^-^i) 
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Of interest is the participation ratio si = Si/S of a given site to the total motion. Its average is 


Si = 1/N. Its second moment, conditioned to the total size S', is easily extracted from (40), 


E{sl\S) = -- 


^{N - l)e^^^^^(cS + ArM;)erfc 


2N^VS 


(41) 


We now stndy the limit of a large number of sites N in Eq. (40). There are (at least) two relevant 
regimes depending on how the driving w scales with N. 


First regime: w = 0(1) (“many avalanches”): Consider the case N ^ oo with w fixed. In this 
case, typical values of S = of order 0{N). Consider S = ' (empirical mean 

avalanche-size Si), which is distributed according to 

N{S-w)^ 


P{S) = exp 

2 


20 FS 


4S 


-^N^oo S{S — w) . 


(42) 


The joint probability P(Si,S), is given by Eq. (40) (with the change of variable S —)■ NS), and 
admits the large-limit 


F(Si,S) ~iv- 


w + cS 

-vexp 


20FSf 

tc(l -I- c) 


1 [w + cS — {1 + c)Si)' 
'4 


3 exp — 
'2 \ 4 


4Si 


1 (1 + c)2 {w - Si)' 


P(S) 


6 {S-w) 


(43) 


2^S( V ^ 

Hence the jump of the center of mass becomes peaked at S = w, while the individual sites keep 
a broader jump distribution. The local avalanche statistics is the same as the one for a particle 
submitted to the parabolic driving force mA{w — Ui), and to the elastic force from the center of 
mass of the interface, c(S — Ui). This observation extends to any number of particles Upart = 0(1) 
with respect to iV: in the large-limit, the particles become independently distributed according 


to the law (43). This picture is the “mean-field” regime usually studied in fully-connected models 


mm, and here derived in a rigorous way. Note that in this case, due to a cancellation in ( [S| ), 
the participation ratio scales as E(s^|S') = 0(1/A^^) which shows that si is typically of order 1/N. 


Second regime: small driving w = 0(1/A^) (“single avalanche”) We now focus on the regime 
w = w/N with w hxed. In this case S = is typically of order 1 and is distributed 

according to 

(S-w)^' 


P(5) = 


w 


2^/TtS''- 


exp 


4^ 


(44) 


We now compute, using (40), the joint PDF of S and Si in the scaling regime Si = 0(1) hxed, 
P{Si,S) 


(45) 


w/N f-il + c)^Si 
-^ exp 




w + cS 


- 3 -exp - - 

2v^(^-^i)i \ 4 


l{w + cS-{l + c){S-Si)y 


4(5 - 5i) 


The hrst factor is reminiscent of the density of avalanches and contains a non-normalizable 
divergence ~ . However (HO) implies a cutoff on small Si of order The scaling w = w/N 
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allows to isolate single (quasi-static) avalanches (in the interpretation of the BFM avalanche process 
as a Levy process discussed above) and the factor of 1/N is the probability that the site i = 1 is 
part of the avalanche. In this regime, the fluctuations are large and the participation ratio scales 
as E{sl\S) = 0{1/N). 

6. Spatial shape in small systems N = 2,3. 

In this section we analyze the PDF of the spatial avalanche shape in the small-driving limit, 
Wi = w ^ 0, mostly for N = 2,3. It already exhibits a saddle-point which allows us to discuss 
the general-iV case below. The analysis can be repeated for finite Wi. Similarities and differences 
give insight into the link between the quasi-static distribution and hnite driving. This is done in 
Appendix F 


N=2,3 We start with N = 2, for which the different models we considered are all equivalent. To 
hx notations, we study the linear chain with PBCs (see Section 2.1) and m = 1. The quasi-static 


PDF of the shape (37), conditioned on the total size S, reads 

2 c 


p{s\S) = 


-c^S 


.2c(1-2«)^ 


s(l-s) 


(46) 


4^/n{s{l - s))2 

We noted s = si = Si/S, the shape variable of the first site. The behavior of this PDF is 
summarized on Figure [Tj For small S, typical avalanches are mainly distributed on one site. As 
S increases, the most probable avalanches become more homogeneously distributed over the two 
sites, and for S larger than Sc = -^, the probability distribution is peaked around s = | and 
the avalanche is extended over the whole system. We call this phenomenon the shape transition: 
For small total size, the most probable avalanches have max(sj) ~ 1, whereas for large avalanches 
max(si) ~ 1/N = 1/2. 


p(s I S) 



Figure 1. Shape transition of the quasi-static PDF (46) for TV = 2 and c = 1 in the linear chain 
with PBCs. For S = 0.15'c (black, solid curve) and S = 0.35'c (blue), the distribution has two 
symmetric maxima. For S = 5Sc, the distribution is peaked around s = 5 (red, upper curve). The 
transition occurs at 5 = S'c = 3/8 (black, dashed curve). 


The case N = 3 for a linear chain with PBC is similar. For S < -^, the quasi-static density 
distribution of the shape p(si, S 2 , S 3 = 1 — Si — S 2 | 5 ') has three symmetric maxima corresponding 
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Figure 2. Shape transition of the quasi-static shape distribution for iV = 3 and c = 1. From left 
to right: S = 0.5; 1; 2. 


to avalanches mainly centered on a given site, whereas for S > ^ there is only one maximum at 
Sj = |. This can be seen on Figure]^ 


General N This study already gives some insight into the structure for generic N: the quasi-static 
distribution of the shape p{^S) exhibits different saddle-points, whose positions and stabilities 
depend on the value of S. For small S, avalanches are preferentially located on a single site j 
and max(sj) ~ 1. As one increases S, the most probable avalanches are more and more extended. 
The analytical calculation of the properties of these saddle points is difficult. However, we can 
generalize the shape transition observed for N = 2,3: The symmetric conhguration dehned by 
\/i, Si = ^ (a situation corresponding to inhnitely extended and uniformly distributed avalanches) 
is always a saddle-point of translationally invariant models. This saddle-point is only stable for 
S > Sc{N), which is computed in Appendix G for the fully connected model, and for the linear 
chain with PBC. The result is 


Sc(N) = 






^N- 


(47) 


;(Ar'5 + 12Ar^ + 0(A^^)). 


160 %^' 

This critical value gives the scaling of the total size above which most probable avalanches are 
uniformly distributed on all the interface. Below this scaling they adopt a more complex structure 
(e.g. they are localized on several sites, possess maxima, etc.). Let us already mention that 
other saddle-points of the shape PDF are numerically studied in Appendix I , where the results 
are compared to the one obtained in the next section for the most probable avalanche shape in a 
continuum model. 
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7. Continuum limit: avalanches of an elastic line and typical shape of avalanches 
with large aspect ratio 

7.1. Avalanche size PDF and density in the continuum limit 

We now study the generalization of the previous result to the continuum Brownian-force model 
with short-ranged elasticity for a line of length L 

'Ijdi'Uxt V Uxt {Ijjxt 'Fxt) “1“ 


Here ^xt is a gaussian white noise with ^xtix't' = d{x — x')6{t — t') and the boundary conditions are 
either free or periodic. Starting from rest at t = 0 and imposing a driving Wxt > 0 for t > 0 such 
that ^^Wxt = Wx, we note the total displacement of the interface Sx = The method used 

in the discrete case can be extended to derive the PDF of avalanches in the continuum. Another 
route is to consider the continuum model as the appropriate N ^ oo limit of the discrete model, 


as is detailed in Appendix H Both procedures give the same result, which, for the dimensionless 
PDF of continuum avalanches, includes a functional determinant 


P[Sx 


1 

2 


det (M) exp 




'0 


''X 

M{x, y) = — ^-i:{V'^)xy + (5(x - y) ( 1 -h 




4^. 

Wx-Sx + ^V^Sx 

Sx 


(49) 


Here is the usual Laplacian, {y^)xy = d''{x — y). Dimensions can be reintroduced as in the 
discrete case using Sm = Sm is the avalanche-size scale of the continuum theory. The hrst 
12 also comes from a determinant and could be included in the dehnition of the 


factor (rr-^ 

1 lx ' 

operator M. 

As in the discrete case, the mean displacement {Sx) satishes —'V^{Sx) + {Sx) = Wx- For 
instance, if the driving is only at one point, Wx = w5{x), one has {Sx) = The case of a 

general Wx is obtained by superposition. This is consistent with the discussion in Section As 
in the discrete case, the mean displacement gives the avalanche shape in the limit of large driving 
(plus an 0{^/w) Gaussian noise). 

One can also study the homogeneous quasi-static limit: w{x) = tc —)■ 0 and S{x) = 0(1) 
uniformly in x. Then P[S] ~ with p[S] the quasi-static density of sizes of continuous 

avalanches, also obtained as the limit of the discrete ones. 


p[Sx 


dxSx)BC[Sx 


exp — 


dx 


(S, - 


(50) 


■ V Jo 4S, 

From now on we set m = 1 (by a rescaling of x). The term BC[Sx] depends on the chosen boundary 


conditions with BC[Sx] = || (resp. BC[S] = for the periodic case (resp. free case). 


Other continuum models Our discrete setting allows us to obtain the avalanche-size PDF of 


various continuous models, Eq. (49) being generalizable to an interface of internal dimension d. 


xy by changing V'^Ux J dyCxyUy. The 


One may also consider an arbitrary elasticity matrix a 
continuum limit of the formula for the PDF of the shape conditioned to the total size, either at 


hnite w, see Eq. (22), or for tc —)■ 0 (quasi-static limit), see Eq. (37), are also easily derived. 
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1.2. Rewriting the probability measure on avalanche sizes 


We now wish to determine the most probable shape of quasi-static avalanches, in the limit L —)■ cx) 
To render the problem well defined, one needs to specify two scales. A natural choice is the total 
size S = f^dx Sx and the spatial avalanche extension (or length) i, i.e. the size of the support of 
Sx- While the avalanche-size PDF P{S) is given by the ABBM result (10), the existence of a finite 
extension i (i.e. local avalanche sizes being strictly zero outside a hnite interval) is non-triviaQ 
Here it naturally arises in the search for saddle-points of the shape PDF: we only found solutions 
which vanish outside of an interval. This property was also shown recently in [15] where the PDF 
of the extension P{1) is computed. 


In the following we study the shape distribution at hxed S and i. We do not take into account 
the term implementing boundary conditions in (50) since it should not play a role in the bulk (this 
hypothesis is explicitly checked on the discrete model in Appendix I). So we write the density of 
continuum avalanches Sx as 



p[S] n dS 

X 

«|S] = 

To eliminate the factor of (n. Sx) in the measure, we set 
= <h2(a;) . 


(51) 

(52) 

(53) 


The integration d<F(x), thus the integral over <F(a;) runs from —cx) to oo. To further 

simplify the calculations, we note that the problem is invariant by translation. We thus impose the 
center of the support to be at a; = 0. This leads to the dehnition of the reduced shape s{x) = (jP^x) 

Sx = ^s{x/t) = ^‘^{x) = ^(j)^{x/t) , [ dx(j)^{x) = l , |a;| > ^ 0(a;) = 0 . (54) 

t t j_i 2 

Note that to study fluctuations around the saddle point it is more convenient to use 0(a;), but the 
saddle point itself can be obtained equivalently using s{x) or 0(a;). Below we use 0(a;), but also 
indicate the corresponding formulas for s(x) when these are simpler. 

We search for the most probable shape in the limit of small driving, at fixed size S and 
extension £. The path integral takes the form 


P[d.^(i')exp (^-^ - Fh4<I>]) 



(55) 


■*" In general the shape of avalanches depends on the driving. However, an avalanche following an arbitrary driving 
(in particular in a quasi-static setting more usual for experiments, see Sec. in the BFM is a sum of quasi-static 
avalanches (Sec. |^, whose spatial structure is, by definition, independent of the driving. 

* In a mathematical sense it may be a peculiarity of the BFM in d = 1 with short range elasticity. Of course 
rapid decay in space is expected more generally beyond some support region of extension £, and often obtained in 
numerical simulations. 
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The boundary conditions are 0(|) = 0(—|) = = 0 and 

J ^ dx(j)‘^{x) = 1 . 

Note the appearance of the factor of ^ in front of the “elastic” energy. 


(56) 


7.3. The saddle point for large aspect ratio 


'he path integral (55) is for large dominated by a saddle-point. To enforce the constraint 


(56), we minimize 'Hei)*/’] —A. J _[/2 with Lagrange multiplier A, leading to the saddle-point 

equations [Jj 


A6(x) = , 50^3;)" _ lOcf^xf^jx) 

^ 2 6(j){x) ^ ^ 0(x)3 0(a;)2 


(57) 


In order to hnd the solution (^o,0o(3^)) of ( |57| ) satisfying the properties written in (54), we hrst 
obtain numerically, using a shooting method, another solution (^i,0i(x)) of (57). We impose 
Ai = 2.5 X 10®, 0i(O) = 1, 0i(O) = 4>i{0) = 0, and look for the correct shooting parameter 0i(O) 
such that the numerical solution has a support of finite size [—Xc, with the desired behavior 
at the boundary, i.e. (j)[{—Xc) = = 0. The obtained (unique) solution has the following 

properties: 4>i{0) = —276.797090676018, Xc = 0.162713, y/0i(a;) ~ 7.85883(a;c — x) for x —>• Xc and 
0^(x)dx = 0.106289. We now take advantage of rescaling, setting 


j2xc 


(58) 


This function is automatically a solution of (57) with a different Lagrange multiplier = (2xc)'^x4i, 

^ ^ ^ -1 

2 i 


and the desired properties (54). By multiplying (57) by (j)o{x) and integrating for x G [—7,1] (using 
())o(±|) = 0), we obtain the relation 77ei[0o] = A 


0 . Numerically we hnd 


So := 77ei[0o] =Ao = (2xe)Mi = 2803.8 ± 0.2 . 


(59) 


An estimate of the numerical accuracy is given. The error is mostly due to the imprecision in 
determining Xc- 

Alternatively, a variational solution can be used. We make the ansatz 


<PvAx) = Wc (x^ - |)M 1 -h ^ Ci(x^ 


i=l 


and 


(60) 


The behavior at the boundary x = ±7 is chosen in agreement with the numerical solution of the 
saddle-point equation. One can also show that this ansatz leads to an energy which remains hnite 
at the boundary. The c-dependent normalization Me is chosen s.t. /_(/2 = 7- ^ given 

vector c = {ci,..., Ci^^^}, one then evaluates ^/[^var]- Using a Monte Carlo algorithm, the minimum 
energy is searched by steepest decent in the space of all c with given = 7- 7^3 

Figure we show that for the shape of the avalanche, this procedure rapidly converges against 


H The saddle point equation has a simpler form in terms of s{x). It reads: ^[s"(x)/s(a;)]" — \ [s'’{x)/s{x)Y = A. 
Hence s”{x)/s{x) is a Weirstrass function which diverges as ~ (x ± Xc)“^ at the boundaries. 
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s(x) 



(5s(x) 



Figure 3. Left: The function so(a::) = as obtained by solving the differential equation (57) 


(red solid curve). This is contrasted to the variational ansatz (60), with one (blue dotted), two 
(green dashed) and 15 variational parameters (black-dashed, indistinguishable from the solution of 
the differential equation). Right: Difference between the solution of the differential equation, and 
the best variational solution. 


the solution obtained by solving the differential equation (57). Our best estimate is for Vax = 15, 
where we find 

c = {-1.00301, 20.6871, 83.4237, 211.353, -270.898,179.973, -72.6636,16.3962, 

-12.2786, 6.11179, -0.33042,11.777, 0.750034, -6.77598, -4.56253} . (61) 

This result is compared to the numerical solution of the saddle point on Figure]^ The energy of 


this solution gives us, in good agreement with Eq. (59), the variational bound 
So < 2803.96 . 


(62) 


In Appendix I we confront this result to a study of the optimal shape in a discrete setting. There 
we also show (see also Figure [Toj below) that this saddle-point is stable. Hence, the reduced shape 

of an avalanche becomes deterministic in the limit of ^ 1: s{x) —)■ so(x) = 0o(^) 


probability one. Formula (55) then shows that Sq is measurable in the tail of the distribution of 
aspect ratios. 


Proba(A/r) 


4^ S/F>1 ( „ S 


exp 




(63) 


with possibly some sub-dominant factors, as e.g. a power-law. This is confronted to numerics 
below. 


7..{. Simulations: Protocol and first results 

Protocol. Here we describe the simulation used to numerically study the shape of avalanches. 
We use a discretization with N = 512 points of the equation of motion for the velocity in the 
BFM (48) using periodic boundary conditions for a system of total size L = N. The mass is 
chosen as m = 10/L in order to get a scale-free statistics for a wide range of events. The other 
parameters are set to unity, r] = a = 1. The time is discretized using a time-step dt = 0.01 and a 
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<Sj > 



Figure 4. Measurement of (Si) and comparison with the exact result (Si) = with 

ig = 256. The total moment is measured as {S) = 99.461. 


discretization scheme identical to [in]. Simulations are done via Matlab and results are analyzed 
using Mathematica. At t = 0 the system is at rest and we choose to drive it using a kick of size 
5w = 100 on a single site. This is motivated by the fact that we want to study (single) quasi¬ 
static avalanches: the value of hw is chosen to be small in adimensioned units —Sw ~ 7.4.10“^. 

cr 

Following the discussion of Section]^ and Appendix D we thus know that an avalanche resulting 
from our driving protocol can either be a “small” avalanche 0{6w‘^) or, with a small probability 
Po = 0{6w) a quasi-static avalanche of total size S = 0(1) (we neglect the 0{6w‘^) probability 
that several quasi-static avalanches have been triggered). Schematically, we write 


PiS) {1 - PoVF'(S) + poPi,{S) , (64) 

where io is the driven site. Here “5” (S) is not a true delta distribution since in the BFM the 
interface always moves, but it rather denotes the PDF of all the small, non quasi-static avalanches, 
which is expected to depend highly on the driving. This is made more precise below, and in 
particular we discuss how we identify the quasi-static avalanches and po from our data set. 

We stop the simulation for the rare events when an avalanche reaches the periodic boundary, 
since we are interested in the distribution of shapes on an inhnite line. For every generated 
avalanche, we numerically compute its shape characteristics S, I (avalanches are indeed observed 
as having a hnite support) and s{x) (discretized with i points). We report results using na = 2.10^ 
simulations of a kick. As a hrst verihcation, we check on Figure a coarse-grained information 
on the spatial structure by measuring the mean local avalanche size. The discrepancy at the 
boundaries can be attributed to the fact that we stop the simulation when an avalanche reaches 
the PBCs. This is the only bias expected in our procedure. It is not a problem since for the rest 
of the article we are interested in observables at large 5/F^, automatically excluding the largest £. 


Consistency check of So = 2804. We predicted above that Sq controls the tail of the distribution 
of aspect-ratios. Numerically, we hnd that this distribution possesses a power-law part coherent 
with an exponent of 2 and an exponential cutoff for large S/C with a prefactor coherent with 
So = 2804: Proba(S'/f'"‘) ~ C/exp{—SoS/C) (see left and center of Figure i- We also remark 
that the exponential cutoff function seems to entirely control the PDF of S/C for “massive” 
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Figure 5. Different histogram of the PDF of S/F obtained numerically with different binning 
procedures for the x axis and scale for the y axis. Left: log-log histogram of the full distribution. 
Center: log histogram of the distribution for aspect ratio S/F > 1/5£q. Right: log histogram 
of the distribution for avalanches of extension £ > 1/m. The black line on the left emphasizes 
the observed power-law behavior VT:o\)a.{S/F) ^ F/S'^. Blue lines are fits using an ansatz of the 
form Proba(S'/£'^) ^ exp{—£oS/F). The red line is a fit using only the cutoff function: 

Proba(5'/.^‘*) ~ exp{—£oS/F). 


avalanches, of extension i > 1/m (see right of Figure |^. Obviously this does not constitute a 
precise measurement of Sq, but rather a verihcation of its non trivial value, which can probably 
only be understood by studying the complete spatial structure of avalanches as we did. 


Identifying quasi-static avalanches. From now on we restrict our numerical results to avalanches 
of extension £ > 10 to obtain a decent spatial resolution. This also allows us to isolate quasi-static 
avalanches. Avalanches with extension larger than 10 only represents 3.5% of the data. Obviously, 
this is not a proof that this subset of avalanches only contains quasi-static avalanches, and one 
needs to check that it has the statistical properties of a set generated by the quasi-static density. 
One “test” is to study the number n^Si of avalanches of total size S larger than Si, for which the 
quasi-static hypothesis implies. 


n>S2 = ^>5i 


LI ds)ds 
L 


(65) 


;s>(S')dS’ 

where p was dehned in ( [26| . Numerically, we hnd that this relation holds for all Si, S 2 larger than 
‘S'min = 0.5 (see Figure [^. We thus further restrict our set of avalanches to avalanches of total 
size S > 5 'min . Note that though our reduced set of avalanches now only contains 2.7% of the 
total number of avalanches, it contributes to 99.44% to the hrst moment (S'). (This gives a precise 
sense to Eq. (64) with po = 0.027). We do not further study the other avalanches here, since their 
characteristics is highly dependent on the chosen driving. 


The convergence to the saddle-point. We now check the striking prediction that the shape of 
avalanches becomes deterministic in the limit of large S/i^. To this aim, we measure the distance 
between the optimal shape so(a^) = 0o(^) the simulated shapes s{x) using either the or the 
(squared) canonical norms (see Figure]^. As expected, we hnd that the mean value of these 
quantities at hxed S/i"^ converge to 0 as S/i'^ becomes larger. However, we hnd that the rate of 
convergence of these quantities is slower than what is expected from perturbation theory (this is 
developed in the next section), which predicts for both a convergence as i‘^/S. This will be taken 
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n(>S) 



Figure 6. Left: n>s measured from the datas (blue dots) and compared to the quasi-static 
prediction ((65), black line) with S 2 ^ S (^i can be chosen anywhere in [0.5,10®] and n^Si is 
measured from the datas). 


(/|s(x)-so(x)|dx > 



{j{s(x)-So{x)fdx ) 



Figure 7. Left: (resp. Right:) Mean-value at fixed S/F of the (resp. squared 
L^) norm between the optimal shape and the simulated shape ^J^^^Ax\s{x) — 50 ( 2 :)| (resp. 

x) — So{x))‘^). Inset: log-log plot of the same quantity, fitted with a power-law {F/S)^F 
(resp. {F/SyFy Error bars are given using a Gaussian estimate and a numerical measurement of 
the variance. The fits with power-laws are of low quality, but sufficient to prove that the convergence 
is slower than F/S. 


into account when comparing the numerical results to the prediction of perturbation theory for 
the fluctuations around the optimal shape. 

The mean shape of avalanches. Finally, we verify on Figure [^that the mean shape (s(x)) is given 
by the optimal shape so{x) for large S/i^. We also explicitly check that the mean-shape decays as 
{x ± 1/2)^ close to the boundaries. The agreement is very good, though one can notice that the 
numerical mean shape is slightly flatter than expected. This observation motivates a study of the 
fluctuations of the shape around the optimal shape. 
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<s{x)> 



<s(x)) 




Figure 8. Left: Mean shape obtained by averaging over the 1000 avalanches with the largest S/F 
(blue dots, 0.0011 < S/F < 0.0041), compared to the optimal shape So{x) (red line). Right: test 
of the predicted behavior s(a;) ~ (a: + 1/2)^ close to the boundaries. 




Figure 9. The coefficients multiplying the different terms in 'H 2 [ 4 >o,S(t'] (left) and 'H 3 [(j)o, 6 (t'] 
(right), after replacing 6 (j){x) —>■ (a;^ — 1/4)^ and 6 (j)'{x) —)■ a;^ — 1/4. This shows that d(j){x) must 
have the same behavior ^ (a:^ — 1/4)^ as (j)o{x) at the boundary x = ±1/2. 


8. Fluctuations around the saddle point 
8.1. Field theoretic analysis 

We now study the fluctuations around the saddle point </>o(x). To this aim, we set 
(j){x) = (j)o{x) + 5(j){x) . 

Expanding the action yields 

^el[0] = ^0 + 2To [ 0o(^)^0(^) +’^2[0O) ^0] + + ••• 


'H2[0O,<^0] = / d(j){x) 


20(j)^^{xy(j)”{x) 15(/)'o(a;) 


41 


± 6 <j)'{x) 


2lO(j)o{xf , .,/// n2 


'H3[0o,^0] = 5 / 6 (j){xyS(j)'{x 


2sixll ^30o(a^)^ - 00 (2:) 00 (2:) 00 (2:) 


0o(a;)2 


+ 77 ^ 0 ' 


00 (a^) 

0o(x)2 


-6(j){x 


(j)o{xy 

3 0'o"(a;)0'o(x) ± 0(;(x)2 


4:6(j){x)6(j)' {x) 


+ sy\xy 

2 0o(^)^ 


00 (a;) 


( 66 ) 

(67) 

( 68 ) 


0o(a;)3 


(69) 
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The first term in Eq. (67) comes from the saddle-point equation (57) at 0 = 0o, = 

\'t>(x)= 4 >o(x) together with (59). We have used our freedom to integrate by part to arrive at 
these expressions: For 7^2 [0O5 we gave a form in which each term is proportional to the square of 
a h0-derivative. For the cubic term, which is used in perturbation theory our strategy is different: 
Since derivatives of {6(f){x)6(f){y))-^^ are numerically unstable, we wrote this expression without a 
second derivative 5(p'\x). 


To evaluate the coefficients, we use the variational ansatz (60), with the optimal c of Eq. (61). 
The plot in Figure shows that (50(a;) should have the same behavior ~ — 1/4)^ as (t)o{x) at 

the boundary x = ±1/2. We therefore make the ansatz 


5(f)[x) = aQVo{x) + ^ \a2n-lVn{x) ± a2nUn{x) 


n=l 


The basis Un{x), Vn{x) is constructed using Gram-Schmidt out of 

2 r 


no(a;) = y - 1 ± cos(27rx) 

Vn{x) = (—1)”'+^ cos(27r(n ± l)x) ± cos(27rx) for n > 1 
in + 1) sin(27ra:) ± (—1)"'+^ sin(27r(n ± l)a;) 


u„\x] = 


f ±n±l 


This basis is orthonormal. In this basis, the energy 7^2 [00, h(f)] can be written as 

1 7 1 X 


(70) 

(71) 

(72) 

(73) 

(74) 


'x,y 




This dehnes M which we now diagonalize. Its lowest eigenvalue is Aq = 2Bo, with eigenfunction 


6 (j)o{x) = (j)oix). This can be proven with the help of the saddle-point equation (57). The higher 


eigenfunctions 5(f)n{x) have n knots, see Figure [loj Since A4. is symmetric they form an orthonormal 
basis. The spectrum is massive (no soft massless modes); we observe that lnA„ ~ 13.1 ± 0.256n, 
i.e. the eigenvalues grow in geometric progression. This ensures that a truncation at rimax = 10 is 
sufficient for practical purposes. 

A delicate problem is to obtain results at hxed = 1. To do so, we write for the 

expectation value of an observable C>[0] 


{o[« = -^ I mown !'f{x) 


exp 




|77el[0O,^0] 




= — j V[cl^]0[cf\ 5 


(a;)-l 


X exp ^-^|772[0 o,50]-7^0 j ± 773[</o, ± 774[0o, ^0] ± •• 


(75) 


We subtracted the constant So from the energy in the path integral and used the constraint 


f^(/>(x)^ = 1 to rewrite the linear term appearing in (67) as a quadratic term: 2So f^(/>o(x)d(/>(x) = 
—So 6(j){xy. It ensures that the minimum of the exponential factor at 6 (j){x) = 0 becomes a 
global saddle point; in addition, the lowest-energy fluctuation 50o has zero energy. If we write 
0(a;) in the basis of eigenmodes 6 (j)n{x) of Ai, i.e. 

OO OO 

(j){x) = (j)o{x) + ^ an6(()n{x) = (1 ± ao)(()o{x) ± ^ an6(j)n{,x) , 


n=0 


n=l 


( 76 ) 
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then 


ct>{xf = 


+ y^^anh(t)n{x) 


n=0 


(1 + ao)‘^ + 


(77) 


n=l 


Solving (/)(x)^ = 1 for oq yields 


ao — 






n=l 


2 tXJ 

ao = -- + ••• 

n=l 


(78) 


With this, the path-integral (75) can be written using equations (76) and (78) as 


{0[(t>]) = -^11 ^ 


n=l 


S r .y ^ A — A 

^ exp I - — { ^ + 773[0O, ^0] + 774[0O, ^0] + •••))• (79) 

n=l 

The factor of (1 — ^n) ^ comes from the derivative of the (5-function, which has been used to 

eliminate the integration over qq. Note that the Jacobian of the transformation from d0(a;) to 
n„da„ is det (J0n(x))^g[_i_i]_„gp^ = 1, since the 6(f)n{x) are orthonormah 

Hence, to leading order in an expansion in the expectation value of an observable of 


6(j){x) can be obtained using the decomposition J0(x) = ^0*(^)’ where ao is given by (78) 

and the a* are centered Gaussian variables with correlation matrix Ai' defined for i, j > 1 by 

t ( 5 , 


{0‘iO‘j)M' 


■ <“> 

One then uses Wick’s theorem for expectation values of (50. As an example, the 2-point 
correlation function is 

°° °° , £8 
(J0(?/)(50(^))^ {aiaj)Mf S4>i{y)S4>j{z) + 0\ 

7=0 j = 0 


S2 


^ 6A{y)6Aiz) ^ 


Xi-Xo 

1=1 


S2 


(81) 


8.2. Generating a random configuration, and exact sampling 

Our setting allows us to generate a random fluctuation with the measure given by the the leading 
behaviour of "H for large S/i^: Denote by gn a series of uncorrelated Gaussian random numbers 
with mean zero and variance 1. Then 


OO 

\ rand ^ ^ 


9r 


=> anJ0n(x), with an = \ -^ ^_ 

n=0 ^ ^ VAn-Ao 


for n > 0 , 


(82) 


and ao given by Eq. ( [T^ . In Figure 11 (left) we show as an example the expectation of (50(x)^ 
(solid blue line). This is compared to the average over 500 realizations drawn with the measure 


(82), repeated 5 times (the three gray-blue lines, lower set of curves). To illustrate the importance 
to properly eliminate the mode 0o(a^), the upper (red) curves are obtained without the constraint 
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6<P„{x) 



Figure 10. Left: The spectrum of A4. The smallest eigenvalue is Aq = 2fo (given with precision 
10“^ for Umax = 10. The next two eigenvalues are Ai = 5.143Ao, and A 2 = 19.20Ao. Eigenvalues for 
large modes grow exponentially with the mode, In A„ ~ 13.1 + 0.256n (black dashed line), showing 
that the spectrum of fluctuations is massive. The lowest modes are colored in red, blue, orange 
and cyan. Right: Plot of the first four eigenfunctions in the same colors as the corresponding 
eigenvalues. S(j)n{x) has n nodes. 


6ip{xf 



Figure 11. Left: Plot of the fluctuations S(j){x)^ (blue solid line), and including the mode 6(j)o (red 
solid line). The dashed lines are averages over 500 samples using Eq. (821, including (top pink) or 
excluding (bottom, blue-gray) this mode. 


on 4>^{x), i.e. including fluctuations proportional to ipoix) (with amplitude ~ I/a/Aq), and not 
constraining them by Eq. (78). 


On Figure 12 we show flve realizations for the shape drawn from the measure (82), and 
compare this to numerical simulations at the same ratio S/F^. The agreement is quite good. 

We can use this formulation for an efficient algorithm, known in the literature as importance 
sampling 120 ]. One writes 


{om 


1 



( 83 ) 














Avalanches in the BFM 


26 


s(x) 



s(x) 



X 


Figure 12. Left: Plot of the normalized shape [(/)o(a;) + \J^5(j){x)\^ for ^ = 3 x 10 Right: 
The same functions from numerical simulations. 


In the second line we reintroduced the full Hamiltonian "Hei using Eq. (67). We will compare to 
simulations below. 


8.3. The leading correction to the shape at large sizes 

For large the mean shape is given by the optimal shape so(a;). For smaller this mean 
shape becomes flatter, an effect which we now investigate using perturbation theory. Consider 

{(5s(a:)) := <s(a:) - M^)^) 

S 


(84) 


x)6(j){x) + 6(j){x) ) ( 1 - ^H 3 [ 0 o, 5(j)] + 


S 


dcpixf) _ -2(j)o{x){6(j){x)'H3[(j)o,d4>] 

/ M' 


M' 


M' 

+ *2(5) ■ 


The notation M! indicates that all expectation values are taken at Sjl^ = 1, making the factors 
of ^ explicit. 


8.J^. Fluctuations of the shape for large avalanches 

We now consider the fluctuations of the shape of an avalanche in perturbation theory: 
((5 s(x) 2>^ := (s(x)2) - (s(x))^ 


[Mxf + 2(j)o{x)6(j){x) + 6(j){xyY { 1 “ ^'^3[0o,<50] + 




M' 


[4>o{xf + 2(j)o{x)d(j){x) + d(j){xy] - ^'H 3 [ 0 o, hep] + ... 

£S2 




M' 




(85) 


Note that the only term which survives is the contraction between one 6(p{x) of each factor s{x) 
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—{6s(x)) 

fA 





Figure 13. Left: (resp. Right:) normalized mean shape displacement {S/F{s{x) — sq{x))) (resp. 
shape fluctuations (S'/£^(s(a:) — So(a:))^)c)- Red line: result of perturbation theory (84) (resp. (85)). 
Dashed-blue line: result from importance sampling using (83) for F/S = 900. Dots: results from 
the simulations for avalanches with aspect-ratio 0.9/1800 < S/F < 1.1/1800 (7023 samples, green), 
0.9/900 < S/F < 1.1/900 (946 samples, blue) and S/F > 1.1/1900 (734 samples, red). We take 
advantage of the symmetry of the observable {s{x)) = {s{—x)) to symmetrize the numerical result. 
We estimate error bars using the difference between the original result and the symmetrized one. 


8.5. Asymmetry of an avalanche 

Another interesting observable is the asymmetry A of an avalanche, dehned by 
A:= 2 X 4>^{x) . 


( 86 ) 


By construction —1 < .4, < 1. The asymmetry has mean zero (4.) = 0, and variance given in 
perturbation theory by 


{A"^) = 16 J xy(j)o{x)(j)o{y) {6(j){x)6(j){y))^, 


r,, = 1.1 X 10-M ^ 


(87) 


8.6. Comparison of the perturbative corrections to the numerics 

We had already shown some results of our numerical simulations above. For large S/C, the 
perturbation theory developed in the preceding section gives the correction {hs{x)) of the mean 
shape to the saddle-point solution, as well as the shape fluctuation ((5s(a;)^)^ around the saddle- 


point. However, as already pointed out in section 7.4, the scaling of these quantities with a factor 


of C/S is not seen in the convergence of the numerical simulations to the saddle point, see Figure 
This indicates that, even at S/C ~ 10“^, the simulations are not yet in the perturbative 
(hrst-order) scaling regime. Non-linear corrections are still important, and ^ {5s{x)) as well as 
p {SsfxY)^ still depend on This is illustrated on Figure 

left of Figure]^, corrections to 


he 


As can be seen on the left of Figure 13 (as well as on t 
the mean shape are very small, of the order of 10“^, difficult to measure, and at the limit of our 


simulations. The red solid line is the perturbative result (84). The points correspond to the same 


quantity from the numerics with increasing S/C from green over blue-gray to red (see caption 
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for the precise parameters). The dashed blue line is obtained for S/t" = 1/900 via importance 
sampling, see equation (Hi El One remarks that the amplitude is lowered as compared to the 
perturbative result, in qualitative agreement with the simulations. In view of the difficulty of 
the numerical simulations, it is very encouraging that at least a qualitative agreement has been 
obtained, and that importance sampling explains why the observed corrections are smaller than 
the perturbative result, in agreement with intuition: the shape has to remain positive. 

The fluctuations around the mean shape, ^ (5s(a;)^)^, are given on the right of Figure [l^ with 
the same color code as previously. One sees that the numerical results approach the perturbative 
result for large S/i'^. In this case, importance sampling predicts fluctuations slightly smaller 
than our numerical simulations, which converge more quickly towards the perturbative result. 
We remark that numerically the estimation of ^ {5s{x)‘^)^ is less sensitive than the estimation of 
^ ((5s(a;)). This may be explained by the fact that only the latter quantity involves non-linearities 
of BL at dominant order in 

For the asymmetry we find ^ = 1.1 x 10“® in perturbation theory, and 5.97 ±0.04 x 10“® 

via exact sampling for S /= 1/900. Numerical simulations give ^ {-A?) = (7 ± 2) x 10“® 
for the largest avalanches F/T^ > 0.002 (37 samples), (5.6 ± 0.3) x 10“® for the data with 
1.1/900 < 5/£^ < 0.002 (697 samples), (4.7T0.2) x 10-^ for the data with 0.9/900 < S/l‘^ < 1.1/900 
(946 samples) and (3.05 ± 0.05) x lO'^ for the data with 0.9/1800 < < 1.1/1800 (7023 

samples). Once again we see that the order of magnitude is correctly predicted (an already non¬ 
trivial achievement), and that the numerical results get closer to the perturbative one as S/^^ 
increases. 

From a conceptual point of view it is interesting to note that most of the amplitude of the 
“double-peak” structure observed on the right of Figure 1^ is due to the first sub-leading mode 
501 (x) with one node at x = 0 (see Fig. 10). The same holds true for {AA). 

In conclusion, we have seen that the numerical results agree very well with the theoretical 
prediction at large and that the mean shape of avalanches is given by the optimal shape 

So (a;) (Figures and [^. The consequence for the tail of the PDF of was successfully 
verified (Figure]^. For finite S'/f"^, namely fluctuations around the optimal shape, we only got 
a partial, though already satisfying agreement: The discrepancy with the perturbative results 
was clearly identified as a consequence of strong non-linearities, even for the largest This 

was qualitatively understood by an implementation of importance sampling, though the remaining 
discrepancy raises the question of wether our simulations are sufficiently precise to measure these 
delicate observables (Figures [I^ and 13). 


9. Application of our results to realistic interfaces and stationary driving 

Up to now we considered avalanches following a stopped driving (see Section [^. However, as 
discussed in [El [H [13] this setting also yields the densities for the statistics of quasi-static 
avalanches in the steady state (Middleton state) for stationary driving in the quasi-static limit 

ft For S/F = 1/900, about 44% of the proposed configurations in the importance sampling have a zero-crossing in 
s{x), and therefore do not contribute. The measured expectation of the weight is (1) = 1.61 ± 0.012, showing that 
averages are not dominated by a few configurations. 
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{wt = V and v —)■ 0^). These are the avalanche densities dehned in Section hence the 
denomination nsed in this article. 

Fnrthermore, it was shown in Ref. HI, that the BFM is the mean-field theory of an avalanche 
in the qnasi-static limit for an interface in short-ranged disorder with eqnation of motion 


^odtUxt / CxyUyt “l- Tn (wxt Uxt) T Uxt) ■ 


( 88 ) 


The disorder-force correlator is given by F{x, u)F{x', u') = Ao(m — u')6^{x — x') with Ao(m) a fast 
decaying fnnction as |m| —)■ oo and Cxy a convex elastic kernel. The prediction of the fnnctional 
renormalization gronp (FRG) for snch systems is that, in the qnasi-static limit, when m —)■ 0 
and for d = due — e, e A 0 (due = 4 for short-ranged elasticity and more generally due = 27 for 
g{q) ~g->.oo the physics becomes nniversal in the small-m limit (e.g. independent of microscopic 
details of the disorder) and entirely controlled by only two relevant conplings, the renormalized 
friction 1]^ and the renormalized disorder cnmnlant Am,{u). The (rescaled and renormalized) 
second cnmnlant of the disorder at the fixed point is non-analytic and exhibits a ensp. It is 
nniformly 0(e), allowing to formnlate a controlled pertnrbative expansion of any observable. For 
observables associated to a single avalanche, it was shown in min] that near the npper critical 
dimension due only the behavior of A^ near zero, i.e. its ensp, Am{u) ~u_j.o= —am\u\ plays a role. 
In this context, the mean-held theory for single-avalanche motion is the BFM stndied here, with 
renormalized parameters 7 —?• 7m and a —)■ exm- Hence, the avalanche densities derived in Section]^ 
are exact for interfaces at their npper critical dimension. They also open the way to a pertnrbative 
calcnlation for d < due- Interestingly, some physical systems described by ( 88 ) are at their npper 


critical dimension, as e.g. domain walls in certain soft magnets for which 7 = 1 |21j . 


10. Conclusion 

In this article we obtained an exact formnla for the joint PDF of the local sizes of avalanches in 
a discrete version of the BFM model. This resnlt is valid for an arbitrary elasticity matrix and 
arbitrary monotonons driving. This allowed ns to derive the densities describing the qnasi-static 
avalanches in the limit of small driving, and to disenss in depth the physical pictnre nnderlying 
this avalanche process. We presented two applications where it was possible to go farther in the 
analytical calcnlation of detailed physical properties. For the fnlly connected model we obtained 
the joint distribntion of the local and global jnmps. This allowed ns to retrieve in a rigorons way 
the nsnal large-A limit, as well as a new regime, and finite-A information. 

We then presented another application by analyzing the most probable shape of avalanches of 
a given size and extension, first for systems made of few conpled particles, then in the continnnm 
limit for an elastic line with short-ranged elasticity. Qnantitative resnlts for the optimal shape and 
the flnctnations aronnd it were obtained and compared to a nnmerical simnlation of the model. 

Let ns conclnde by stressing that, since onr formnla was obtained in a general setting and 
contains all the spatial statistics of avalanches, it shonld be possible to extract from it a variety 
of new information on their spatial strnctnre of direct experimental interest. It wonld also be 
interesting to compare onr resnlts for the shape of avalanches to other models throngh simnlations 
or experiments, the BFM being the relevant mean-held theory for varions more realistic systems. 
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Appendix A. Recall of the result for the generating function 


For completeness, we recall in this section, the derivation, here in a discrete setting, of the exact 
result for the generating function of the BFM ([^. Related derivations can be found in [131 E]- 
The original equation of motion, including the quenched noise term dtFi{uit) reads 

N 

r]dtUit = ^ Cijiijt - - Wu) + dtFi{uit) . (A.l) 

j=i 

We use the dynamical field theory formalism [221 [23] which allows to compute the disorder average 
of any physical observable 0[ii\. We introduce N response fields Uu such that disorder averages 
can be computed as 


0[u\ = / 


(A.2) 


The dynamical action splits into a deterministic, quadratic part and a disorder part: S[u, u] = 
S'o['U,'u] + Sdis [«,«], with 


5*0 [ft, m] = 


N / N \ 

Uit[ rjdtUit - Cijiijt + m^{uit - wu) 

i=i V j=i J 

j^Uit CijUjt + 


N 

E 

2=1 


m UitWit 


(A.3) 


where in the second line, we made an integration by part assuming ii vanishes at infinity. The 
disorder part of the action is 


5'dis [ft, U\ 





(A.4) 


it contains all the correlation of the Gaussian force (|^. As noted in [T31IT1] . the action functional 
can be simplihed using the Middleton property recalled in the main text, valid for our setting: 
t2 ^ ti < > Uit2 Y Uit^ SO that 


dtdt'\uit - Uit'l = Uitdt'sgn{t - t') = -2uit6{t - t') . 


(A.5) 


This leads to 


Sdis [ft) "tt] 


= —a 


N 

y~! / ■ 

i=l 


(A.6) 


It is straightforward to check that the replacement dtFi^uu) —)■ y/2auit^l used in the main text 
leads to the same action. This shows that both theories are equivalent for this choice of initial 
conditions. As written, the action is linear in ft: this simplifies the calculation of the generating 
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functional of the velocity field G[A,w] = e^i=i 


N 


2 r — 


JJ 5 I Aa + aul + T]dtUit + ^ CijUjt 
it \ j=l 


m^Uit 




(A.7) 


In the last line, the response held is solution to the “instanton” equation [lansiiii] 


N 


Xu + au^n + TjdtUit + ^ CijUjt - m^Uit = 0 . (A.8) 

j=i 

It is imposed by the delta functional. Note that this evaluation involves a w-independent Jacobian, 
which equals unity since we have supposed the interface to be at rest and stable for t < 0, so that 
if Wit = 0 then iiu = 0. The above result is thus correctly normalized. Equation (A.8) must 
in general be supplemented by some boundary conditions, depending on the observable (e.g. if 
Xu = 0 for all i and t > ti, we should also have itu = 0 for all i and t > ti). Note that a rigorous 
version (in discrete-time, without path integral) of this result was given in [13]. In the main text 
we are looking for the statistics of avalanches Si, which is obtained using constant sources Xu = A,, 


and for which one can look for constant solutions Uu = Ui of (A.8). 


Appendix B. Tests of the main formula, computation of moments and numerical 
checks. 


We checked (18) using two methods: the hrst one consists in solving exactly the instanton equation 


for small values of N in an expansion in powers of c for a given elasticity matrix. This gives 
an approximation of the Laplace transform, which can be inverted to give the joint probability 
distribution up to a certain order in c. This program has been successfully achieved up to O(c^) 
for N = 2, O(c^) for N = 3 and O(c^) for N = 4. The other method consists in numerically 
computing various moments of the probability distribution, which can then be compared to the 


exact results that use the instanton equation (|12|): the cumulants are given by 

d 




dXi 


^1 


-lnG(A)) 

/ Ai =0 


N 


k=l 


dvk 


and theses derivatives are numerically computed using = J-X where Jj,- = —2vi6ij -|- C, 


dXi, ■ ■ ■ dXu 


T-l 


(B.l) 


Vi=0 


dXj 






as 


seen from (12). 


Appendix C. Backward Kolmogorov method for a kick driving 


In this section, we provide another verihcation that (18) is correct when the system is driven by a 


kick (i.e. Wu = Wi6{t)). For simplicity, we directly consider the dimensionless equation of motion 


N 


dthit ^ ^ ^ij^jt ^it T Wu T x/ 
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N 


-t>0 




(C.l) 


j=i 


where in the second line we used the definition of Cij 


and wrote the equation for f > 0 when 
Wit = 0. For a kick, it is equivalent to consider the equation of motion with Uit=o = 0, or to consider 


the equation without driving for f > 0 (C.l) supplemented with the initial condition Ui^t=o+ = un- 


The generating function G is still given by G(X) = F or a kick, we can write it as a 

conditional expectation value on the process without driving (C.l): G(X) = G(X, w, 0, oo) where 
G is defined as 

G(A, w, ti, tf)=E( 




(C.2) 


where iiu evolves according to (C.l) for all times and E |^... = Wij denotes the average on the 

stochastic process without driving ( C.l| ) conditioned to the initial condition = Wi. We now 
derive a partial differential equation (PDF) fo G, similar to a Backward Kolmogorov equation, 
using a splitting of [ti,tf] into + ht] U [ti + St,tf] with 6t small: 


G{X, w, ti, tf) = E dtuu+j:Zi Ai 4*+'* dtu,. 


^iti 




Uiu = tc* + o{6t) 


(C.3) 


Where in (C.3) we used that uu is continuous. The expectation value in (C.3) can now be split in 


two parts. We can first average over the noise for t G [ti, ti + St], with St small, or equivalently on 
the velocity variation Swi := iii^ti+st — Ui,ti = Ui^ti+st — Wi, as obtained from the equation of motion 


(C.l). Secondly, we average over the noise in [U + St,tf] (these are independent) knowing that the 
velocity at t^ + St is Ui^a+st = Wi + Swi, i.e. 

G{X,w,U,tf) = E|,^q ( E ( 




’Mu = Wi e^»=i + o{St) 


= E{5^„j (^G(A, w + Sw, ti + Sti,tf) iiiti = Wi^ + o(St) 


(C.4) 


The average over {Swi} can be computed at first order in St using Ito’s lemma (we use Swa = 


St 'X2f=i ^ajWj and Sw^^ = 2waSt + 0{St‘^)). This leads to 


N 


^ ^ ^ 1 \ f)G 

G{X,w,ti,tf) = + —(2u;„) ) + St— 

0 = 1 \ ^ j = l ^ 


' dti 


(C.5) 


N 


X 


1 + ^ XiStwt I + o{St). 


2 = 1 


We also expanded the last term at first order in St. In the r.h.s. of (C.5), all generating functions 
are taken at the same position G{X, w, ti,tf). Now the l.h.s. is of order 0{St^) and in the l.h.s., we 
exactly computed the 0{St) term. This shows that the generating function G solves the following 
PDE: 


dG 

dti 


N 

E 

0 = 1 


dG 

dWrv 


N 


d^G 


^ ^ ^ajWj T Qiu‘2 Wa T X(yWi 


i=i 


70 I ^' 0*^0 


(C.6) 
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which is also equal to as a consequence of the time translation invariance of the Brownian 

motion. The initial condition is G{X,w,ti,ti) = 1. 

To study avalanche sizes, we consider the long-time behavior of G to obtain G = G{X, w, ti, oo). 
In this case we can assume that G reached the stationary state, i.e. 


N 

E 

OL=l 


N 


^G r-i d^G . 

y GajWj + XaWa 1 = 0 


dw, 


i=i 


dwl 


(C.7) 


This is automatically satisfied if G is given by (|^ and if the Ui satisfy the instanton equation ([^. 
This provides a connection between the two methods. 

An interesting feature of this method is that one can now write a PDE directly for the 
probability distribution P{w, S) of avalanche sizes in the BFM model following arbitrary (positive) 
kicks Wit = Wi6(t). This equation reads: 


N 

E 

a=l 


N 


dP d‘^P 

y ^ CajWj T^Wn 


dw. 


i=i 


dwl 


Wn 


dP 


= 0 . 


We need to find a solution which satisfies the following boundary condition: 

N 

P{w = 0,S) = l[6{Si) . 


(C.8) 


(C.9) 


2=1 


Let us now discuss its solution. Inspired by our result (18), we make the change of variable 
P{w, S) = F{x, S) with X = w — G ■ S. The equation for F then takes a very simple form: 


N 

E 

a=l 


w. 


d‘^F dF 


dxl dSa 


= 0 


(C.IO) 


where Wa = Xa + GajSj and we used that C is a symmetric matrix. In this new variables, 
we write our main result (18) using the following decomposition: 


F(f , S) = det {Mij)F{x, S) 

, N 


Mij — Gij + 6ij — 




exp 


2=1 



(C.ll) 

(C.12) 


This decomposition sheds some light on the structure of ([I^ , here rewritten as F in (C.ll): it is 
simple to see that F defined in (C.12) already solves (C.IO), F can indeed be interpreted as the 
PDF of the position Xi at ’’time” Si of N independent particles diffusing from the origin at time 
Si = 0. However the result F = F would not satisfy the boundary conditions (C.9). We now check 
that the extra factor det(M) provides the proper solution. In order for (C.ll) to also solve (C.IO), 
the determinant must verify 


N 

E 

a.=l 


Wn 


d‘^ det{M) ~ ddet{M) dF ddet{M) 
~F + 2-- 7 :; -- F 


dxl 


dXrv dXr 


dS^ 


= 0 . 


(C.13) 


Using ^ = —f^F, this implies an equation for det(M) 


2S, 

N 

E 

OL=l 


w. 


cI^det(M) Xad det (M) 9det(M)\ 

dxl Sa dXa dSa J 


(C.14) 
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The first term ^ is equal to 0, since Xa only appears in the a-th column of M. The remaining 

terms vanishes since M depends on Xa and Sa only through the combination This completes 


the proof that our result (18) indeed solves the PDE (C.8). The boundary condition is now satished 
since is a continuous PDF on positive variables and we know (see Section]^ and Appendix 

B) that {Si) = vanishes when Wi —0. 

Appendix D. Poisson-Levy process for normalizable jump densities 


Center of mass We already discussed in the main text the inhnite divisibility property (28) of 
Fw(5'). Given this property, one would like to interpret an avalanche S as the sum of n iid 
elementary avalanches Si with n drawn from a Poisson distribution and s* drawn from a given 
distribution (this dehnes a Poisson-Levy jump process, see e.g.[I7j). This interpretation is valid 


at the level of the moments of Pw{S) (see (30)) but we now show that it does not extend to the 


probability itself. Let us hrst assume that the jump density p appearing in (30) is normalizable 


(see also the discussion in Appendix J). Then one can write p{s) = pop{s) with p a regular 
function normalized to unity J dsp{s) = 1 and po fhe density of avalanches; i.e. the mean number 
of quasi-static avalanches occuring in response to the total driving w is Pqvj. Using the following 
identity: 


dsi ■ ■ ■ ds„(e^^i - 1) ■ ■ ■ - l)p(si) ■ ■ ■ p(s„) 


m=0 


(Pqw)" (-Pqw)^-'- 
m! [n — m)\ 


dsi ■ ■ ■ ■ --Pism) 


(D.l) 


(30) can be rewritten as (performing the sum over n > m): 

j dSe^^P^{S) = f2^-^^e-P°'^j dsi---ds^e"(^i+-+^-)p(si)---p(sJ . (D.2) 


m=0 


This leads to a formula for the probability, P^iS) = X]m=o e~^°'"(p*)’^(5'). Here (p*)^ 

denotes m convolutions of p with itself, making the interpretation in terms of a Poisson jump 
process transparent. One can dehne the “complete” avalanche-size density as 

dP^{S) 


p{S) = 


dvj 


.=o = -po6{S)+p{S) . 


(D.3) 


Where here the hrst equality holds in the sense of distributions. This total density appears 
as the sum of the regular density p{S) (dehned in the main text) and of a delta singularity 
that accounts for the hnite probability that the interface does not jump. As a consequence, 
|w=o = = J dS'e^'^p(S') = J dS (^e^^ — l) p{S). For the ABBM model, the scale 

invariance of the Brownian motion leads to an accumulation of small avalanches of arbitrary small 
sizes, leading to po = cxd (in particular for any w > 0, Pw{S = 0) = —)■ 0) and one can not 

dehne p. The formula |w=o = f dS (e^^ — l) p(S) is however still valid and allowed us to 


prove (30). 
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Levy Process for the interface The generalization to the interface is immediate: in this case, the 
LT of Pw{S) reads 


/ OO 

^_ r\ /j 


Wi^... Wi^ 




(D.4) 

lb'. 

n=0 

where the second sum is for all (ii,... ,in) € A^}"' and the Vi variables are functions of A 


solutions of (12). Using our conjecture (^33L we obtain 


/ OO 

^—n fi 


Wi^... Wi„ 


n=0 


n\ 


nd"fi,(e“‘‘ - l)ft.(4) ■ ■ ■Pi„(4) 


(D,5) 


i=l 


which is the multidimensional generalization of (|30j) and shows that the densities Pj{S) entirely 
control the moments of Puj(S'). It is also in agreement with the interpretation of an avalanche S 
as a superposition of independent avalanches, as already discussed in the main text. 

Appendix E. Details on the fully connected model 

Here we detail the calculations leading to the results of Section and give some results for the 
fully-connected model driven by a single site. 

Marginals distributions for uniform driving For uniform driving, the matrix C and M entering 


in (18) admit the following simple expressions, allowing us to evaluate det M in a concise way: 


Cij = (1 + c)5ij - ^ , Mij = 5ijj{w + cS/N) - ^ 


N 


1 


N 


det M = w{w + cS/Nf-^ JJ — , S = 


2=1 


(E.l) 


fc=l 


This leads to (|38|). Various marginals of this PDF can be computed by noting that the Laplace 

(E.2) 


transform of entering into (38) reads 


dsp^,s/N{s)e-P^ = 


We write the joint PDF of local and total size as 


N 


p(ys) = i p(s). 


(E.3) 


2 = 1 


For any 1 < m < A — 1, the marginal Pd^i,..., Sm}, S) can be computed as 


w 


N 


p{{Si,...,s„,},s) =—^^——Y\p^^s/N{Si) / TT Pw,s/N{Si) 

W + Cb/Jy-^-^ /v-Af q. 

2=1 2_-^i=m+l i=ri 


w 


w + cS/N J-1 
' 2=1 




N—m)w,{N—m 


2=m+l 
m 

,s/A-(S- 

2=1 


(E,4) 


Where the multiple convolution of Pw,s/n{.s) has been easily calculated as a consequence of the 


simple structure of it’s Laplace transform. In particular, this leads to the formula (40) of the main 
text. 
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Single-site driving Taking Wi to be non-uniform breaks the permutation invariance i t-)- j of the 
problem, making the computation more complicated than for the uniform case. Another solvable 
case is tCj = 0 for i ^ 1, for which the PDF ( [I^ takes the form 

N 


P{S) = 


Siwi 


S{wi + cS/N) 


Pwi 


,s/n{Si) 


(E.5) 


i=2 


The computation of marginals involving an integration over some Sj for j > 1 is identical to the 
uniform driving case and leads, for l<m<A^ — 1, to 

Q m m 

Po,S/7v(5'j)po,(V-m)S/v(*5 - Y1 

In particular, we obtain 


(E.6) 






20F(5 - ^i)3/2 


exp 


(m;i + cS/N - (1 c)^i)' 


X exp 


In this case S = Yld=i typically of order 1 and is distributed according to 


((iV-l)(cF/iV)-(l + c)(^-Fi)r 
4(5 - 5i) 


45i 


9{S-Si). 


P{S) = 


Wi 


2y/nS'- 


exp 


AS 


(E,7) 

(E.8) 


The large-A^ limit now exhibits a single non-trivial regime, with wi = 0{N^), and for which (E.7) 
admits the limit 

c5i ( {wi - (1 + c)5i)^ 


F(5i,5) = 


Wi 


2^sl 


exp 


X exp 


(c5-(l + c)(5-5i))' 


45i 


9{S-Si). 


4(5-5i) 

Remarkably, in this case one can even integrate over the total size to hnd the marginal PDF P{Si 
in the large-A^ limit. 


PiS,) = 


Wi 


exp 


{wi - (1 + c)Siy 


(E.IO) 


2^5f V 

In agreement with the physical intuition, this is the ABBM result for a particle with driving 
rn^{wi — u) and c(5 — u), as discussed above, and 5 = 0, since the center of mass has not moved 
appreciably. 


Appendix F. Shape for small N at finite driving 

Here we briefly discuss what becomes of the shape transition observed in the quasi-static PDF of 
avalanche shape at hxed total size 5 of the linear chain with PBCs (see Section when one is 


interested in the full PDF for hnite Wi = w as given in (22). For N = 2 and w < there is now 
an additional regime with two transitions instead of one: 
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• S < x/3V3 16CW+3 . distribution of s is peaked around 

• - 8 g-»^-v^V 3 -i 6 e '»>+3 ^ g ^ - 8 c^+v^V 3 -i 6 ct^+ 3 . distribution possesses two symmetric maxima 
around s = |. 

• S > -8c^^+v^V3-i6cw+3 ^ retrieves a single maximum at s = |. 

The first regime is new, and was not captured by the study of p. For small tc —)■ 0 it corresponds 
to avalanches smaller than the lower-scale cutoff S < which are not described by p as we 
know from Section In this regime, the fact that the saddle-point again corresponds to uniform 
avalanches with s = 1/2 is not a consequence of elasticity (as noted in Section]^ local avalanche 
sizes are even independent in this limit), but is related to the fast decay of po{s) at its lower cutoff 
(see Section 1^. For larger tc > the intermediate regime disappears, and the most probable 
avalanches are homogeneously distributed. Indeed, as w increases, the motion of the interface 
becomes mostly deterministic and the remaining fluctuations become negligible. 

The case = 3 is identical. For tc < ^ the finite w probability distribution exhibits the same 
three different regimes with boundaries 0 , ^^d interpretation is 

identical to iV = 2. 


Appendix G. Stability of infinite, uniform avalanches. 


In this appendix, we compute the value Sc{N) such that avalanches uniformly distributed over all 
the system, and of total size S > Sc{N) are stable. We do this for the fully-connected model and 
for the linear chain with PBC s, for which uniform avalanches uniformly distributed are always an 
extremum of the quasi-static density p (for uniform driving fi = 1 ). As such, Sc{N) is the value 
of S above which all the eigenvalues of the hessian of the quasi-static distribution at this uniform 
saddle-point are negative. Since this saddle-point and the elasticity matrix are translationally 
invariant, the Hessian of the logarithm of the probability at the saddle point is a circular matrix 
given by 




log p{^S) 


S 

Ys' 


i=s — + h, 


1 

2 ^ 2 ' 


a/3- 


(G.l) 


dSadSfS 

c is the elasticity matrix of the model (here = 1), s = 1/A is the uniform local avalanche size at 
the saddle-point and ha/s depends on the chosen model as h^js = —jj^ + j^{4:6a/3 + Sa,i 3 -i + Sa,i 3 +i) 
for the linear chain with periodic boundary conditions, and haB = + jihaB for the fully 

connected one. The eigenvalues of these matrices can be computed using a discrete Fourier 
transformation, showing that they are indexed by a wave-vector q = ^ with fc = 1,..., A — 1. The 
g = fc = 0 mode does not intervene since it corresponds to a uniform displacement of the interface, 
which is forbidden by the fact that we work at hxed S\ dsi = 0. The eigenvalues of the Hessian 
are all identical for the fully-connected model: A/.c. = ^ 2 ^ hor the linear model they 

are given by \q = —^[1 — cos(g)]^ + ^ + 2cos(g)]. In the latter case, the most unstable 

mode is g = leading to the following critical values 


S^N) 

^pbc(a) 


3A 

9 ’ 


A ,1 1 , ^ /2vr,,, 

2c2(l-cos(f))A2 + 


(G.2) 
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16c^7r^ 


+ 12N^ + 0{N^)) 


(G.3) 


Appendix H. Continuum limit 

Here we detail the scaling that allows to find the probability distribntion of the dimensionless 
continnnm avalanches P[S' 3 ;] knowing the probability distribntion of the discrete case P{S). We 
denote for clarity the continnnm field as Ut{x), x G [0,L], and its iV-point discretization as 
Uit = Ut{i^). We will add indices c and d to distinguish between physical quantities of the 
continuum and discrete models. An easy way to ensure that the statistic of the discrete case 
corresponds to the statistic of the continuum one is to compare the different terms in the dynamical 


action (see Appendix A) : 


The disorder term: dx acMxfMx) ^ Eti 


N 


The elastic 

-N L 


term: 

L 


It 

J2ZiIt^itCd{ui+it - 2,Uit + Ui_u) = dx j^Ut{x)ccA.Ut{x) ~ 






• The driving term: Eii It mluuWit = fj' dx I mlut{x)wt{x) ~ Eii It 

This indicates that the quantity of the discrete model should be q = ^Cc 

and ad = In particular, the rescaled quantities which appear in the text, in the formula 

for the dimensionless discrete distributions are -% = and St, = Note that we will 

rrij mi m L m 

choose everywhere in the main text q = 1. This implies that the probability distribution of the 
dimensionless rescaled continuum avalanches denoted by Pc is given in terms of its discrete analog 


P = Pd given in (18) as (introducing the explicit dependence in the driving): 


F,|S(x),u.(x)] = Jim ^ - 


N 




(H.l) 


where here S = (S'(Li/A))j=i^...^ 7 v and w = {w{Li/N))i=i^...^N. This leads to the formula of the 
main text. Note also that for //-dependent observables, one should choose rjd = ^rjc- 

Appendix I. Optimal shape in the discrete model 


Here we compare the results on the continuum optimal shape with the discrete case. This is not 
only a consistency check, but also allows us to compare the results of the optimization when we 
include boundary conditions, and to investigate the stability of the shape. We choose to work 
on the discrete model with an elastic coefficient set to unity, which corresponds to a A-point 
approximation of the continuum model with a line of length L = N, i.e. the index i of the discrete 
model is the coordinate of the continuum line (see Appendix H). In the continuum, the optimal 


reduced shape sq is obtained for total size S and extension i fixed, and contains all the probability 
when 3> 1. To compare this result with the discrete model we used two different optimization 
procedures on the discrete probability. We always impose the total size S and optimize on the 
shape variables Sj = St/S with 

(i) either the two central points tuned to coincide with the optimal continuum result: we note 
Umid the integer part of A/2 and impose = |so(0.5/£). 
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Figure II. Comparison between the most probable shape of length £ = 10 with iV = 12 computed 
using optimization on H (blue dots) or p (red dots), using procedure (i), and for different total sizes 
S from left to right: ^ = 10“^, 10“^, 1. The influence of boundary conditions quickly decreases as 
S/F is increased. 


(ii) either N — I successive shape variables fixed to be small (below we use s* * = 10 ® ) 

Procedure (i) is an indirect way to impose the extension by imposing that the avalanche shape is 
peaked around some region, whereas procedure (ii) is closer to the continuum setting where we 
directly imposed the finite extension. In both cases we impose S' 3> to obtain a true maximum. 
The optimal shape is always found to be symmetric, which allows us to impose this condition to 
study reasonably large N. The result of the optimization is then compared with the prediction 
from the continuum theory: Si = =s>£4 jSo{i/i). One can then 

• Verify that the optimization on p (including boundary conditions) or Ft alone (defined in the 


continuum in (51)) give the same results. It is already obvious for i N and Figure II 


explicitly shows that it is always true for S' even if £ ~ V. This validate the hypothesis 

made in the continuum that boundary conditions do not play a role for large S/i'^. 

Using an optimization on FL, we can verify that the discrete optimal shape coincides with 
the continuum one. The results are shown in Figure [1^ One can see that, apart from some 


discretization artefacts, procedure (ii) give results in agreement with the continuum result. 
On the other hand, procedure (i) leads to a shape with an effectively larger extension. This is 
in agreement with the idea that the property that avalanches have a strictly finite extension 
is only a feature of the continuum limit, as explained in Section 7^2, and is coherent with the 
idea that procedure (i) only imposes a “characteristic” extension in the discrete setting. 

Finally, we can study the behavior of the maximum eigenvalue Xmax of the Hessian of the 
discrete Hamiltonian FL at the most probable shape (since the eigenvalues are negative it is 
the maximum one that is the closest to 0 and that controls the stability of the saddle-point) 
using procedure (i). The behavior of the eigenvalues of the Hessian with S is trivial: since 
S can be factorized in front of the Hamiltonian, they are proportional to S. However, in 
the discrete case, there is no way to see the scaling ^ emerge from the Hamiltonian. Still, 
we clearly numerically find (see Figure 13) that Xmax scales with for £ —)■ 0. This thus 
provides an alternative verification that the saddle-point is stable, and that it’s stability is 
controlled by S'/£^ 3> 1. 
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Figure 12. Most probable shape in the discrete model obtained using numerical optimization on 
H with procedure (i) (blue dots) or procedure (ii) (red square) with = 30 and £ = 16 (left) or 
£ = 22 (right), compared to the continuum saddle-point prediction so{xl£)ji (straight line). 
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Figure 13. Maximum eigenvalue of the hessian of the hamiltonian at the numerical optimum as a 
function of -p for large, fixed S with procedure (i). 
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